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ABSTRACT 



High levels of exozodiacal dust are observed around a growing number of main sequence stars. The origin of such 
dust is not clear, given that it has a short lifetime against both collisions and radiative forces. Even a collisional 
cascade with km-sized parent bodies, as suggested to explain outer debris discs, cannot survive sufficiently long. 
In this work we investigate whether the observed exozodiacal dust could originate from an outer planetesimal 
belt. We investigate the scattering processes in stable planetary systems in order to determine whether sufficient 
material could be scattered inwards in order to retain the exozodiacal dust at its currently observed levels. We 
use N-body simulations to investigate the efficiency of this scattering and its dependence on the architecture of 
the planetary system. The results of these simulations can be used to assess the ability of hypothetical chains of 
planets to produce exozodi in observed systems. We find that for older (> lOOMyr) stars with exozodiacal dust, 
a massive, large radii (> 20AU) outer belt and a chain of tightly packed, low-mass planets would be required in 
order to retain the dust at its currently observed levels. This brings into question how many, if any, real systems 
possess such a contrived architecture and are therefore capable of scattering at sufficiently high rates to retain 
exozodi dust on long timescales. 

Key words. Planetary systems, Comets: general, Planets and satellites: dynamical evolution and stability, Planets 
and satellites: general, Zodiacal dust, Infrared: planetary systems, Methods: numerical 
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1. Introduction 

Excess emission in the far-infrared or sub-mm, asso- 
ciated with cold, dusty debris discs is regularly ob- 
served aroun d hundreds of main sequence stars (see 
IWvattl d2008l) for a review). There are now a grow- 
ing number of stars with similar emission, but in 
the mid-infrared from material at higher tempera- 
tures. These observations have been associated with 
dust very close to the star, reminiscent of the so- 
lar system's exozodiacal cloud. There are now sev- 
eral dozen such sources discovered either using near or 
mid-infrared, interferometric observations with, for ex - 
ample, FLUOR/C HARA(e.g. lAbsil et al.l l2006l l2008b. 
IONI C/IOTA(e.g iDefrere et all \2QUh, MIDI/VISIR 
(e.g. ISmithetal1l2009bh. or using their mid-infrared 
spectra with Spitzer doiofsson et al.l 1201 2t iLisse et all 
I20T2L [2009h . Such sources have been compared to the 
solar system's exozodiacal cloud, but in general are sev- 
eral orders of magnitude brighter, for example 77 Corvi 
is 1250 ± 260 times brighter than the solar system's 
exozo diacal dust cloud at 10;um dMillan-Gabet et aU 
.These observations can be modelled to deter- 
mine estimates on the geometry, mass, position and 



grain properties. A well studied example is the 450Myr 
old A -star Vega (lAbsil et al.H2006l |2008t [Defrere et al.1 
1201 ll) . where the best fit model found using Bayesian 
statistics to assess a grid of models, suggests that the 
emission is dominated by small grains (< 1/um) in a 
10' 9 M ffi ring between ~ 0.1 - 0.3AU dPefrere et al l 
1201 ll) . More detailed compositional information can be 
determined for the handful of brighter sources detected 
by their Spitzer spectra i n the mid-infrared, for exam- 
ple lOlofsson etail (120121) : ILisse et al.1 (1201 2l) . amongst 
others. 

Once the potential alternative of a stellar compan- 
ion is ruled out for the interferometric observations, 
the only explanation for this excess emission is emis- 
sion from high levels of small dust grains. Such small 
grains, however, have a short lifetime, both against col- 
lisional processes (years) and radiative forces (years). 
If the observed discs are to survive for longer than year 
timescales, the small grains must be replenished. In 
outer debris discs, the observed small grains are thought 
to be continuously replenished in a steady-state colli- 
sional cascade from a population of unseen, large par- 
ent bodies. However, mod e ls for s uch collisional evo- 
lution (e.g. IWvatt & Dentl (|2002|) ; IWvatt et al l (1 19991 
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2007bb: IStark & Kuchnej J2009I) ; iKrivov et alj d2008b ; 
Lohne et alJ ( 2008 ) amongst others), find that close to 



the star, the planetesimal belt required to supply the ob- 
served levels of small dust must be so dense that even 
kilometre sized parent bodies are rapidly destroyed by 
collisions dWvatt et al.ll2007ah . Thus, it is not possible 
for the observed small grains to be replenished by an in- 
situ steady-state collisional cascade on long timescales 
(hundreds of years at 1AU). The suggestion has been 
made that the dust could be resupplied in a steady state 
collisional cascade if the populatio n of bodies were o n 
highly eccentric (e > 0.99) orbits dWvattet all 1201 Oh . 
The origin of such a configuration is, however, difficult 
to conceive. Therefore, in general another explanation 
for the origin of the observed small dust is required, in 
particular for the older systems with exozodi. 

Many suggestions have been made in the liter- 
ature to explain these exozodi. The young systems 
(~ 10 - 20Myr) with hot dust are easier to ex- 
plain and it seems probable that we are witnessing 
increase d collision rates during terrestrial planet for- 
matio n dSmith et al.ll2009al l2012t iKenvon & Bromlevl 
120041) . whereas for the older systems, the observa- 
tions are more puzzling. An obvious origin for the 
material is further out in the planetary sy stem where 
it can survive on longer timescales ("e.g. [ Smith et alJ 
12009b!: iLisse et al.ll2012t lAbsil et al.ll2006h . With the 
growing number of detections of planets and/or plan- 
etesimal belts, it seems that outer planetary systems 
are common. In fact some of the observed exo- 
zodis also have detections for out e r planetesimal belt s 



(e.g. Fom alhaut lAbsil et all d2009h:lAcke et a 
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amongst many others). It remains, however, to be de- 
termined how this material is transported inwards and 
whether we are observing a transient or steady-state 
phenomena. 

Recent dynamical modelling of our solar system's 
exozodiacal dust cloud found good evidence that much 
of the observed dust originates from Jupiter Family 
Comets, scattered inwards from the Kuiper belt. Once 
inside Jupiter's orbit the dust's evolution is dom- 
inated by radiative f orces, Poynting-Robertson-drag 
dNesvornv et alJ 120101) . Both scattering and PR-drag 
are transport processes that could occur in extra-solar 
planetary systems. Poynting-Robertson-drag is most in- 
fluential in low luminosity systems, for example, it 
has been suggested to mo del observations of warm 
dust around Epsilon Eridani dReidemeister et ai]|201 ll) . 
Scattering can occur both in a stable planetary system, 
such as the Jupiter Family Comets in our solar sys- 
tem, and at increased rates after a dynamical instability, 
such as the solar system's Late Heavy Bombardment. 
Evidence that we may be observing the aftermath of 
a collision between two large bodies exists from the 
presence of impact processed sili ca in the spectra of the 
hot dust dLisse et al. l l20Lll2009h . Models of collisions 
between large bodies, such as the earth-moon forming 
collision, show that hig h levels of dust are produced 
dJackson & Wvattll2012l) . Such collisions are by nature 
rare, in particular their probability decreases with the 



system age and thus, their probability must be com- 
pared to the growing numbers of systems with exozodi- 
acal dust, in particular, th e number of old (gre ater than 
Gyr) systems (e.g. t Ceti ddi Folco et al.ll2007l) ). 

In this work we investigate whether sufficient mate- 
rial can be scattered inwards from an outer planetesimal 
belt, in a stable planetary system, to resupply an exo- 
zodi, or whether it is necessary to invoke further expla- 
nations. We consider the simple scenario where there 
are (several) planet(s) interior to an outer planetesimal 
belt and use N-body simulations to determine how ef- 
ficiently they scatter material inwards from the outer 
belt. To assess the feasibility of this scenario, these ef- 
ficiencies must be compared with observations of inner 
and outer belts. We start by discussing the details of our 
simulations in Sec. [2] Secf3]presents the illustrative ex- 
ample of single planet systems, which is expanded upon 
in Sec |4] with simulations to determine the scattering ef- 
ficiencies in multi-planet systems. The implications of 
our results and their limitations are then discussed in 
Sec. |5] before conclusions are formed in Sec. [6] 



2. Our simulations 

The main purpose of this work is to investigate the rates 
at which material can be scattered from an outer plan- 
etesimal belt into the inner planetary system, with the 
view that it could contribute to the observed exozo- 
diacal dust. This is necessarily dependant on the ar- 
chitecture of the planetary system; that is the orbits 
and masses of the planets, and the radial location and 
masses of the planetesimal belt and exozodi. We aim 
to constrain the planetary system architecture required 
in order to scatter material inwards at sufficient rates to 
produce the observed exozodiacal dust. We anticipate 
a high dependence of the scattering rates on the sys- 
tem architecture, thus, we hypothesis that the diversity 
of observed systems could potentially be explained by 
a diversity of planetary system architectures. For sim- 
plicity, we limit ourselves in this work to the scenario 
of an outer planetesimal belt. 

In order to assess the validity of this hypothesis, 
we use N-body simulations to determine the efficiency 
of scattering in planetary systems of different architec- 
tures. Given the huge diversity of planetary systems 
found by current observation^], it is necessary to limit 
our parameter exploration to a tiny region of the avail- 
able parameter space, which is chosen such that the 
simulations can be used to make further deductions re- 
garding the remaining parameter space. Our simula- 
tions determine the rate at which material is scattered 
inwards as a function of the total mass of material in 
the outer belt, that is the total mass in all bodies be- 
tween small yt/m grains and km-sized bodies. The mass 
of the outer planetesimal belt can be determined, al- 
beit in a model dependent manner, from observations 
of outer debris discs. For systems where an outer belt 
is not detected, this merely places an upper limit on the 
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mass that could be hidden in an outer belt in this sys- 
tem. Using these masses it is possible to assess whether 
this hypothesis can be applied to real planetary systems 
with both outer planetesimal belts and exozodi. This re- 
quires a knowledge of the rate at which material needs 
to be supplied in order to sustain the observed levels of 
exozodiacal dust. 

Observations of exozodiacal dust can be used, again 
in a model dependent manner, to determine a mass for 
the observed dust. It is clear that the observed dust is so 
small (/im) that it has a short lifetime, either against col- 
lisions or radiative forces, although the exact lifetime is 
difficult to calculate. Multiple epoch observations, sep- 
arated by longer than th e lifetime of the small dust, of 
a handful of systems (e.glCiardi et al.ll200lt lAbsil et all 
120061 120081: IDefrere et alJl201 lh .do not find significant 
variation in the observed fluxes. Therefore, it seems rea- 
sonable to make the assumption in this work that the 
small dust should be retained at its currently observed 
rates. It is, however, not clear whether this small dust 
is resupplied directly, or via collisions between some- 
what larger (maybe only mm sized) grains also orbiting 
in the disc. In order to avoid a complicated discussion 
of the collisional evolution of material in the exozodi, 
we decide to make the base-line assumption that at least 
the currently observed mass of dust must be supplied to 
the exozodi, on the lifetime of this small dust. In gen- 
eral, this is likely to be (significantly) lower than the 
rate at which material with the size distribution of the 
outer belt (where the mass is dominated by large km 
sized bodies) needs to be scattered inwards, since it ig- 
nores the inefficient conversion of this material to small 
dust (discussed further in Sec. [5}. However, if stable 
outer planetary systems cannot even scatter material in- 
wards at this rate, then it is highly unlikely that they are 
capable of replenishing the exozodi. As our fiducial ex- 
ample we take the Vega system, where a best fit model 
finds 10~ 9 M e in between 0. 2yum and 1mm grain s, with 
an estimated lifetime of lyr (IDefrere et al.ll2.01lh . Thus, 
to retain Vega at its currently observed levels, the exo- 
zodi must be resupplied at a rate of at least 10~ 9 M ffi /yr, 
although uncertainties in these values should be taken 
into account in the assessment of our simulations. 



2.1. The scattering process 

In this work our focus is on the fate of planetesimals 
that are scattered during the planet's sculpting of the 
planetesimal belt. Of particular interest are those parti- 
cles that are scattered into the inner planetary systems 
at late times, as these have the potential to replenish an 
exozodi. 

The scattering of planetesimals by planets is a sub- 



erature (e.g. Ravmond et al. 
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20001 Beust & Morbidelli 
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Bolevetal. 20121 


Bonsor & Wvatt 2012; Debes etalj 2012). Planets are 



often invo ked to sculpt the inner edges of debris 
discs (e.g. lOuillen & Faberi 120061: IChiang et al.l 120091: 
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Fig. 2. The time at which particles are scattered inside 
of 1 AU as a function of their initial semi-major axis, a,-, 
for the same simulation as shown in Fig.Q] No particles 
that started outside of a, = 25a c / MOJ are scattered. 

disc of planetesimals, as shown in Fig. [1] will scatter 
particles that are close to it. These particles lie in a zone 
directly surrounding the planet, whose size can be esti- 
mated analytically by consid ering the overla p of mean 
motion resonances, given by dWisdomlll980t) : 



Sa c haos = Cfl 



2/7, 



pi, 



(1) 



iBonsor et al.l 1201 11) . A planet placed at the edge of a 



where fi — the planet mass and C is a constant, 
with various values between C — 1.3 dWisdomlligS 1 
C = 1.5 dDuncan et al~l ll989t lOuiilen & Faberj|20Qi 
to C = 2 dChiang et al.ll2009h in the literature. In this 
work, we use the value C = 2, as we find this is a good 
fit to our simulations (see Fig.[TJ. Almost all particles in 
the planet's chaotic zone will be scattered, given suffi- 
ciently long timescales, however, some particles outside 

of this zone will also be scattered by t he planet. 

Al though simulations such as lOuillen & Faberi 
d2006l) show that planets generally clear their chaotic 
zone on less than 10,000 orbital periods, the effects 
of this clearing can continue on significantly longer 
timescales. Even in our own solar system, after Gyrs of 
evolution, there is evidence of Kuiper belt objects scat- 
tered by Neptune, for example, Jupiter Family Comets 
or Neptune's scattered disc. The best way to understand 
some of these processes is to consider a very simple ex- 
ample. 

In Fig. Q] we present the results of a simulation in 
which a IM^ep planet was placed on a circular orbit at 
2AU, with an initially 'cold' (e < 0.1 and / < 10°) plan- 
etesimal belt, extending from 2 to 3AU. Particles that 
start inside of the chaotic zone, and close to its edge, are 
scattered by the planet onto orbits of higher eccentric- 
ity, in general at either constant pericentre or apocen- 
tre. The scattering timescale is proportional to distance 
from the planet, thus the chaotic zone is cleared from 
inside outwards. In the left-hand panel, it is clear that 
after 10,000yrs of simulation, the chaotic zone is still 
being cleared; all particles initially in the inner regions 
have been scattered, whilst there are still some parti- 
cles on unperturbed orbits in the outer regions of the 
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Fig. 1. The semi-major axes and eccentricities of particles scattered by a single Neptune mass planet on a circular 
orbit at 2AU around a 1M Q star, after 10,000yrs (left-hand panel) and lMyr (right-hand panel). The set-up of 
the simulation is described in Sec. 12.11 Trojan's of the planet are labelled by red crosses. After 10,000yrs, only 
particles within the inner regions of the chaotic zone (Eq. Q]with C = 2, shown by the dashed lines) have been 
scattered, whilst by lMyr, the full zone has been cleared. Particles are scattered onto orbits at constant pericentre 
(q) or constant apocentre (Q), thus lie between the dotted lines q ~ a p i and q ~ a p i - 5a c i mos or Q ~ a p i and 

Q ~ a p i + SUchaos (Q ~ dpi). 



chaotic zone. On the other hand, as shown in the right- 
hand panel, after lMyr of simulation, the full chaotic 
zone has been cleared. Particles are scattered onto or- 
bits at constant pericentre or apocentre, such that the 
scattered particles lie between the lines q ~ a p i and 
q ~ a p i - 5a chaos or Q ~ a pt and Q ~ a pl + 5a chaos , 
where they are scattered multiple times, until eventu- 
ally they are either ejected, or hit the planet or star. This 
can be on significantly longer timescales than the clear- 
ing of the chaotic zone. Such particles are equivalent to 
Neptune's scattered disc in our solar system. 

In terms of the formation or replenishment of ex- 
ozodi, we are interested in particles that are scattered 
close to the star, inside of a limiting radius that we call 
f Z odi- Whether this refers to the location of the exozodi 
or the radius at which Poynting-Robertson or other drag 
forces dominate is irrelevant to the current simulations. 
Particles that are scattered inside of this radius origi- 
nated in (or near to) the planet's chaotic zone, but have 
been scattered multiple times by the planet, generally 
evolving at constant apocentre, until their eccentricity 
is sufficiently high that they are scattered inside of r ZO£ /;. 
Thus, the total amount of material scattered generally 
depends on the size of the chaotic zone and critically 
the planet's mass. 

Fig. |2] shows the time at which particles are scat- 
tered inside of r zo< u = 1AU, by the same IM^ep planet 
on a circular orbit at 2AU, as a function of their initial 
position. The broad trend is for an increase in the time at 
which particles reach r :oc n = 1AU with increasing ini- 
tial semi-major axis. This was anticipated, as the planet 
scatters particles close to it more quickly. This trend 
is, however, very smeared out, in a way that would not 
be the case if we merely plotted the time at which the 
particles were initially scattered. This is because once 
particles have been scattered a few times by the planet, 
their further scattering looses its dependence on the par- 



ticle's initial orbit. Given that most particles are scat- 
tered many times before they end up inside of r ZO di, the 
trend in Fig.|2]is somewhat blurred. 



Critically, Fig. [2] shows that the scattering rates, 
even at late times, can be influenced by particles that 
started on orbits very close to the star. This is unfor- 
tunate since it is not clear whether or not particles are 
expected on such orbits after planet formation. As dis- 
cussed in further detail in Sec. 13.31 the end conditions 
of planet formation are not well constrained. It is, how- 
ever, clear that planets clear a zone directly surround- 
ing them, although there are no precise constraints on 
the size of this zone, which may well differ between 
terrestrial planet and gas giant formation. It is not the 
purpose of this work to investigate the end conditions 
of planet formation, nor to rely greatly upon their out- 
come, therefore, we restrict the initial semi-major axes 
of our test particles to «„„■„, where «„„•„ > a p \. This 
removes particles from very close to the planet, that 
would have almost definitely been accreted onto the 
planet, but ensures that we include particles close the 
edge of the planet's chaotic zone. Given that our initial 
test simulations show that such particles are only scat- 
tered on timescales greater than 10, 0007^/, it seems un- 
likely that all of these particles would have been signifi- 
cantly influenced, by the smaller planetary core, to have 
been accreted during planet formation. Thus, although 
the value of a mm may be somewhat arbitrarily chosen, 
it should nonetheless provide a good representation of 
the behaviour anticipated in a real planetary system. 
This, however, ignores the effect of migration during 
planet formation, that could potentially leave gaps in 
the planetary system. These affects and the influence of 
our choice of «„„■„ are discussed further in Sec. 13. 31 
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Fig. 3. A diagram to illustrate the set-up of our simu- 
lations. The large circle represents the planet, at a p i. 
The real belt, at late times, is shown by the thick line 
between r be it,m = a p \ m + 5a chaos to r be[tm , = ^hi. 
To simulate this belt we need only consider test par- 



lpl,out 



+ 



and 



tides that start between a mi „ 
CLmax - 0-pi,out +?-5a c i mo!S , represented by the dotted, thick 
line. 



2.2. Details of the simulations 

Our simulations are run usin g the RMVS3 integrator a s 
part of the SWIFT package ilLevison & Duncan|[l994l) . 
This has an adaptive timestep. We consider planetary 
systems that have an outer planetesimal belt, running 
from rbeit.m t° rbeit,out, and a total mass of Mbeit in all 
bodies (between the blow-out size, Di,i and a maximum 
size, D„ mx , as in Eq. [5]). The planets are on circular, 
coplanar orbits, interior to this outer belt. We consider 
that the outer planet, orbiting with a semi-major axis of 
c-plout has sculpted the inner edge of the outer belt (see 
for example Fig . [TJ . The outer planet is therefore placed 
interior to the outer belt's inner edge by its chaotic zone 
I.e. Clpl^out — f belt, in - 6a c haos- For our initial sim- 
ulations we only consider a single planet, whilst later in 
this work (Sec. [4} we consider chains of several equal 
mass planets, of mass M p i, separated eq ually in units of 
their mutual Hill's radii, giv en by (e.g. iRavmond et all 
120091: IChambers et al.ll 19961) 



Rh = 



(ai + fib) ( Mi + M 2 
3M» 



1/3 



(2) 



where a\ and fib, M\ and Mo are the semi-major axes 
and masses of the two planets. 

Our outer planetesimal belt is composed of N test 
particles, whose initial semi-major axes are randomly 
selected between «„„„ and a max , initial eccentricities be- 
tween and e„ mx = 0.1, initial inclinations between 
and I max = 10°, and whose initial longitude of as- 
cending node, argument of pericentre and mean anoma- 
lies are between and 2n. Although we consider our 
simulations representative of a belt between rt, e it,m and 
rbeit.out, it is not necessary to simulate particles covering 
the full belt (/.e. a mi „ = r he i,j„ anda max = r belUout ). Initial 
test simulations (see Fig. |2| found that a single planet 
has little to no influence on particles outside of twice 
the width of its chaotic zone. Thus, to reduce comput- 
ing time, we fix a max at a P i out + 26a c i wos . As will be 
discussed in further detail in Sec. 13.31 the planet will 
already have cleared a zone surrounding it during its 
formation. The exact size of this zone is not well con- 



strained and particles this close to the planet are scat- 
tered on average on shorter timescales, therefore we ar- 
bitrarily fix a m i n at a P i out + , although the depen- 
dence of our simulations on this value should be noted. 

The results of our simulations are scaled to the full 
belt of width Tbdtin to r belt, out- Care must be taken when 
referring to the mass of the outer belt. In this work 
we refer to the total mass of the outer belt as the to- 
tal mass in a belt between fbeit.in — Q-pi.out + &@chaos 



and The 



T<r b eh,i, 



, except for in the examples where 



the outer belt radial extent is known from observations. 
This outer radius for the belt is an entirely arbitrary 
choice, and our results can be readily scaled to any 
value. It, however, enables a comparison between ex- 
amples with different outer planet masses to be readily 
made. 

The details of the set-up used in our simulations are 
shown in Fig. [3] and all the important parameters are 
listed in Table 12.21 In order to avoid unnecessary com- 
puting, we consider that particles have left the planetary 
system if their eccentricity increases above 1 or they are 
scattered outside of r max , which we take as 10AU in the 
one planet simulations and 1000AU in the multi-planet 
simulations. 

To determine the rate at which material is scattered 
inwards, we track the rate at which particles are scat- 
tered inside of r zol a, the nominal location of the exozodi, 
which we generally take to be 1 AU. Given the necessity 
to limit the number of test particles in our simulations, 
we consider that each test particle represents a fraction 
of the mass in the outer belt, contained in a suitable 
number of particles that reflect the size distribution of 
the outer belt. The fraction of the outer belt mass rep- 
resented by each test particle is assigned based on the 
particles initial semi-major axis, assuming a mass sur- 
face density of the belt given by l.(r)dr oc r~ a dr, where 
or = 1, and that a ~ r (valid for the low initial parti- 
cle eccentricities used here). Given the potential differ- 
ence between the simulated belt and the actual belt (see 
Fig. [3]), each particle is assigned a fraction of the total 
mass in the belt, extending between rbeitjn and rb e it,out- 



•■part 



1 



r belt,out r beltA 



PP I \ &max ^mm 



(3) 



where N pp is the number of test particles placed be- 
tween a„,j„ and a„ mx at t = 0. Using this formulation, 
this rate can readily be scaled to an outer belt of a differ- 
ent width (i.e. with a different rb e it, ut- The rate at which 
material is scattered inwards is then calculated from the 
fraction of the belt mass assigned to each particle di- 
vided by the time between adjacent collisions. 



2.3. The survival timescale of massive outer belts 

As can be logically expected, and as will be quantita- 
tively investigated in our study, the scattering of bodies 
from an outer belt to an exozodi is facilitated for mas- 
sive outer belts. There is, however, a theoretical limit 
as to how massive such a belt can be, imposed by its 
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Table 1. A list of the parameters used in this work 



Parameter 


Nominal value in this work 


Description 


A 




The separation of the planets, in units of their mutual Hill's radii 


3&chaos 




The radial extent of the chaotic zone in units of the planet's semi-major axis 


p 


ZlOOkgin 1 


The density of the individual particles in the outer belt 






The initial semi-major axis of a particle, before it is scattered by a planet 


a P i 




The planet's semi-major axis 


a p!,oul 


' belt, in 3&chaos 


The outer planet's semi-major axis 


a 


„ , 5a chaos 
pi .out 


The minimum initial semi-major axis of test particles in this simulations 


a'mlx 


^ pi .out ^^^chnos 


The maximum initial semi-major axis of test particles in this simulations 


a fin 


5AU 


Planets are placed between this radius and the outer belt. 




60km 


The size of the largest body in the outer disc 


dr 


7=2 


The width of the outer belt 




d.i 


The maximum initial eccentricity of test particles in this simulations 


ej 




The eccentricity of a particle before it is scattered by the planet 






The disc flux 






The stellar flux 


(e) 


0.05 


The average eccentricity of particles in the outer belt 




10° 


The maximum initial inclination of test particles in this simulations 


M, 


1M G unless otherwise specified 


the stellar mass 


Mdel 




The minimum mass of a planetesimal belt such that it can be detected 






The maximum mass of a collisionally evolving planetesimal belt 


Mpart 




The fraction of the total outer belt mass assigned to each test particle 


M pl 




The planet mass 


N PP 




The number of particles used in the simulation 


Qo 


1507/tg" 1 


The dispersal threshold 


(jfmin 




The minimum pericentre of a particle, with a Tisserand parameter of T,, 
scattered by a single planet on a circular orbit 


Rh 




Hill's radii 


r 




The central radius of the outer belt 


r bell, in 


Qpl,out "t" &&chaas 


The inner radius of the outer belt. 


K belt, out 




The outer radius of the outer belt. 




1AU 


The position of the exozodi, or the inner radial limit within which most 
scattered material will go on to produce an exozodi 


K max 


10 or 100 AU 


The distance outside of which particles are considered removed from the simulations 


' min 




The minimum radius of a planetesimal belt in order that a given mass of material can survive 
against collisions for the age of the system 


1~belt,min 




The minimum radius required by a planetesimal belt such that it can scatter material inwards at 
a sufficiently high rate to retain the observed levels of exozodiacal dust 


T, 




The initial value of the Tisserand parameter 



collisional evolution. This concept of a maximum pos- 
sible belt mass for a given system's age is a very impor- 
tant constraint when discussing our results. It is there- 
fore important to introduce it before presenting and dis- 
cussing our simulations. 

Debris discs are collisional system in which larger 
bodies are continuously collisionally ground down to 
smaller bodies and the smallest bodies are removed by 
radiation pressure. Thus, over time the total mass in 
the disc is depleted. Care must be taken in referring 
to the total disc mass, which here refers to the total 
mass of all bodies in the collisional cascade, extend- 
ing from the observed small dust grains (yum) to km 
sized bodies. There may be larger bodies that are not 
yet in collisional equilibrium, it is however, impossi- 
ble to obtain any information about such bodies from 
the observations and we therefore choose to disregard 
them in this work. Models for the colli s ional evolutio n 
of such discs from IWvatt etal] d!999l) : IWvattl d2003l) : 



IWvatt et al] d2007allbh find that there exists a maximum 
total disc mass that can survive against collisions on 
long timescales. This maxim um mass as a funct ion of 
the disc radius, r, is given by (IWvatt et al.ll2007al) : 

M max (t) = 5.2 X lfr 13 M ffl X (4) 
(p/kgm- 3 ) (r/AU) 13 ' 3 (f XeyJkg- 1 ) 5 / 6 (D ma Jkm) 
(MJMv)- 4 ' 3 ) (e)- 5 ' 3 (f/Myr)- 1 , 

where p is the bulk density of the bodies, r the cen- 
tral radius of the belt of width dr, Q* D the dispersal 
threshold, D„ mx the size of the largest body in the disc, 
M, the stellar mass, (e) is the average eccentricity of 
particles in the disc, assumed to be constant and t is 
the age of the system, or strictly the time after the disc 
was stirred. Written in this form, it was assumed that 
the disc contains a Dohnanyi size distribut ion of bodies 
(iDohnanvil 19691) with a power index of ^ dTanaka et al.l 
[19961) " between the smallest size body that is not blown 
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Fig. 4. The minimum central radius of a planetesimal 
belt, such that a total mass of either M c n sc = 1CT 5 , 1CT 3 , 
1, or 10M ffi , can survive, as a function of age. This is 
calculated using Eq. [6] assuming a steady-state colli- 
sional equilibrium, with bodies up to a maximum size 
Dmax = 60km and fiducial values for the other param- 
eters taken from fits to observations of main sequence 
debris discs, quoted in Sec. 12.31 The results scale di- 
rectly with and the time axis can be scaled by a 3. One planet 



disc mass to the maximum diameter particle is fixed, 
rather than the exact disc mass. Planetesimal belts in 
steady state can exist above and to the left of the mini- 
mum lines plotted. 

Fig. |4] confirms our earlier assertion that exozodi 
cannot be supplied by a steady-state collisional cas- 
cade of massive, unseen, parent bodies orbiting at their 
radial location. Consider the example of Vega, where 
the best fit to the observations of the exozodi find 
1.9 x 10~ 9 M m of mater ial, in 0.2//m to 1mm sized bod- 
ies dDefrereetal.ll201ll) . If this material were to be sup- 
plied in a steady-state size collisional cascade from un- 
seen bodies extending up to D ma x = 60km in size, at 
the same radial location (i.e. inside of 1AU), an abso- 
lute minimum on the total disc mass required would be 
10~ 5 M ffl . Fig. |4] clearly shows that such a mass of ma- 
terial cannot survive at 1AU for longer than ~ 30Myr, 
less than 10% the age of Vega. However, such a mass 
(and s ignificantly more ) can survive for 455Myr age of 
Vega dYoon et alJ2010h PL in a belt outside of60 AU (the 
observed debris disc radius (iMiiller et alJl2010h \ 



factor involving the stellar mass, f^ 4 3 - 



out of the system by radiation pressure and the signifi- 
cantly larger maximum size for bodies in the disc, D max . 
The mass is dominated by the largest bodies. It also as- 
sumes that particles can be catastrophically destroyed 
by particles significantly smaller than themselves, that 
the disc is thin (of width dr) and does not spread radi- 
ally. 

This maximum mass is a strong function of the disc 
radius. We can use it to determine the minimum radius 
of a planetesimal belt in which a particular mass can 
survive against collisions, up to a given age. Thus, the 
minimum radius at which a disc of mass, M,;, JC , can sur- 
vive is 

r min = 6.8 x 10 2 AU x 

(M disc IM @ fl l \M t IM e ) A l l \e) 5 l l \tlMyrfl n 

t /kmr 3,a {p/kgm- 3 r 3,a (^r m (Qr D /Jkg- 1 )^' 26 



The simplest possible planetary system is one contain- 
ing an outer planetesimal belt and a single interior 
planet on a circular orbit. Although possibly not a likely 
structure for many real planetary systems, this provides 
a simple to understand, illustrative example of the scat- 
tering processes that occur in more complicated multi- 
planet systems, which we defer to later in this work. 
We investigate the rate at which particles are scattered 
to the position of the exozodi. However, a single planet 
is not capable of scattering particles infinitely far and 
therefore it is useful to firstly consider analytically the 
position of the planet relative to the outer belt and the 
exozodi. 



(A, 



This minimum radius is particularly interesting be- 
cause it enables us to place a radial limit on the posi- 
tion of an outer belt if sufficient mass is to survive for a 
given age. This minimum radius is plotted in Fig. [4] as 
a function of time. Formally this is the time since the 
collisional velocities in the disc reached high enough 
values to be catastrophic and the system was stirred, 
which we assume to be approximately equal to the age 
of the system. We take fiducial values for the model 
parameters that come from fitting observa tions of de- 
bris di scs around main sequence A stars (IWvatt et alJ 
l2007bh . Thus, the maximum body is D max = 60km, ec- 
centricity (e) = 0.05, the density p = 2700kgm 3 and the 
width — = 0.5. In order to allow some of these parame- 
ters to be easily varied, the time is plotted multiplied by 
a factor involving the stellar mass and the ratio of the 



3.1. An analytical limit on the maximum semi-major 
axis of the planet 

If particles scattered by a planet are to produce an exo- 
zodi, the planet must be capable of scattering them suf- 
ficiently far in, to r zot u. If we consider planets on cir- 
cular orbits, the conservation of the Tisserand param- 
eter places an an alytical limit on how fa r a planet can 
scatter a particle dBonsor & Wvattl [20121) . This can be 
expressed as a minimum allo wed pericentre for the par - 
ticle's orbit, given by (Eq.3 of iBonsor & Wvattl (I2012h 1 



-Tf + 2Ti + 4 - 4 ^JT^Ti 



T 2 - 



(6) 



where T, is the initial value of the Tisserand parame- 
ter, for a particle with semi-m ajor axis, q,-, eccentricit y, 
e, and inclination, I, (Eq.3 of dBonsor & Wvattll2012l) V 
This limit is only strictly valid for low mass planets on 
circular orbits, where long range interactions are rare. 



2 although note the factor of ( 2.1M e ) 4/3 necessa ry in Fig. [4] 
to account for the mass of Vega dYoon et alj|2010h 
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It, however, provides a good approximation to the be- 
haviour even of Jupiter mass planets on mildly eccentric 
orbits, although there may be a small flux of particles 
scattered inside of this minimum pericentre. 

This minimum pericentre can be used to determine 
whether or not a particular configuration of planet and 
planetesimal belt can scatter any material to the posi- 
tion of the exozodi, r ;oi /,. For a single planet it pro- 
vides a limit on the maximum semi-major axis of the 
planet. This depends on the initial orbital parameters 
(Tisserand parameter value) of the particles before they 
are scattered by the planet and is determined by equat- 
ing q„,m(Ti) = r zo di- This maximum semi-major axis is 
plotted in Fig.|5]as a function of the particle's Tisserand 
parameter value, except that the particle's initial eccen- 
tricity is used as a proxy for its value of the Tisserand 
parameter. This means assuming values for the parti- 
cle's initial inclination of 7, ~ 2e,, and initial semi- 
major axis of a, - a p i + 5a c i wos {\MN ep ), which places 
the particle initially at the edge of the planet's chaotic 
zone. In this plot, we consider a Neptune mass planet, 
although the dependence on the planet mass of the re- 
sults is negligible. In terms of our single planet simula- 
tions, given that we start particles in a fairly cold belt, 
with e max = 0.1, Fig. [5] shows that the planet must be 
inside of a p i < 2.2AU if it is to scatter such particles in 
as far as an exozodi at r zoc n = 1AU. 

This analysis can also be applied to the inner planet 
of a chain of multiple planets, as considered later in this 
work (Sec. |4|. In this case, the particle's inclinations 
and eccentricity will have been excited by the exterior 
planets before it reaches the inner planet. For example 
all particles with e, ~ 0.2 by the time they reach the 
inner planet, can be scattered as far as 1 AU by a planet 
at 5AU. 



3.2. Dependence on the planet mass 

Given the analytical limit on the ability of a planet 
outside of 2.2AU on a circular, planar orbit to scatter 
particles inside of r zoc n = 1AU, we consider the fidu- 
cial example of different mass planets on circular or- 
bits at a ; ,/=2AU. The simulations are set-up in the man- 
ner described in Sec. 12.21 and the rate at which mate- 



rial is scattered inside of r zoc n = = 1AU is tracked, 
as a function of the total mass of a belt spanning be- 
tween r MtJn = a p i + 6a chaos and r be it,out = = 3AU. 
Each cross on the plot in Fig.[6]represents an individual 
test particle in our simulations, at the time at which it 
reached r zo di. 

An analytical understanding of the dependence of 
scattering process on planet mass can be found by con- 
sidering the variance in the size of the chaotic zone 
(Eq.[Tji with planet mass. Thus, an increase in the total 
amount of material scattered with planet mass is antici- 
pated. Generic arguments tell us that higher mass plan- 
ets scatter on shorter timescales. Thus, a significantly 
increased rate of scattering, at least at early times, is an- 
ticipated. This can be seen in Fig. [6] Although the high 
mass planets are good at scattering particles inwards at 




0.5 



Fig. 5. An analytical limit on the maximum semi-major 
axis of a planet, if it is to be able to scatter particles 
as far inwards as the position of the exozodi at r zol n, as 
a function of the particle's initial eccentricity. This is 
calculated using Eq. [6] and the particle's initial eccen- 
tricity is a proxy for the Tisserand parameter value, as 
described in Sec. 13.11 This limit only strictly exists for 
planets on circular orbits and vanishes for e, > 0.46 (or 
Tj < 2), or for very low eccentricity particles that do 
not cross the planet's orbit (e, > 0.04 or T, < 3). 
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Fig. 6. Comparing the rate at which particles are scat- 
tered inwards by Neptune, Saturn and Jupiter-mass 
planets on circular orbits at 2AU. The details of the 
simulations can be found in Table [2] Higher mass 
planets scatter more material inwards, but lower mass 
planets continue to scatter material inwards on longer 
timescales. 



high rates on short timescales, on longer timescales, low 
mass planets are better suited to scattering particles in- 
wards at lower rates. 



3.3. Dependence on the initial conditions of the 
simulation 

The precise outcome of N-body simulations depends on 
the exact initial conditions used. In our case, the initial 
orbital parameters of the planetesimals in the outer belt 
can affect the ability of planets to scatter them inwards. 
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Fig. 7. The change in the scattering rates with the width 
of the simulated belt, i.e. a mm and a max , for a lM^ ep 
planet at 2AU, as described in Table [2] The crosses 
represent individual scattering events, whilst the dotted 
lines are straight line fits to the data. There is a negligi- 
ble difference between scattering rates when a max is in- 
creased (green and black points), however, a decrease in 
a m j„ causes a slight increase in the scattering rates, even 
at late times (red or blue crosses compared to black or 
green). 
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Fig. 8. The scattering rates of particles that started with 
different initial eccentricities and inclinations. The sim- 
ulations are all for a single Neptune mass planet at 2AU 
and are detailed in Table [2] The crosses represent in- 
dividual scattering events, whilst the dotted lines are 
straight line fits to the data. There is an increase in scat- 
tering rates with the particle's initial eccentricity or in- 
clination. 



Ideally the initial conditions for our simulations would 
be those found at the end of planet formation process. 
It is, however, not fully clear what these are. There are 
many uncertainties in our understanding of the planet 
formation process and its outcomes, and it is not the 
purpose of this work to investigate these. We therefore 
choose to consider an idealised scenario in which the 
planets are fully formed at the start of our simulations, 
and the planetesimal belt is situated outside of their or- 
bits. 



As planets form they accrete material from their di- 
rect surroundings. Thus, a deficit of planetesimals very 
close to the planet would be expected at the end of 
planet formation. However, we do not anticipate that the 
planet fully clears its chaotic zone, given that our initial 
simulations (Fig. [TJ showed that some particles in the 
outer regions of the chaotic zone were only scattered af- 
ter 10,0007^/, i.e. longer than expected timescales for 
planet formation. As discussed in Sec. 12.21 we therefore 
did not include any particles between a p i and a mm = 
a p i + l5a ^'"°-' in our simulations. Fig. Q shows the the 
range of change that could be anticipated in our scat- 
tering rates, if we have either under or over estimated 
a m j„. At the late times of interest, even the maximum 
possible change is less than an order of magnitude. 

In order to reduce the necessary computation time, 
we also limited the initial semi-major axis range of test 
particles to less than a max = a p i + 28a c i wos . Fig.|7]shows 
that this does not alter the results more than the statisti- 
cal variation between individual simulations. For multi- 
planet systems, there is the potential that three-body 
resonances or secular terms may affect particles much 
further out in the disc, however, test simulations found 
this not to be significant (Sec. [5}. 

Arguably the most significant initial condition for 
our simulations is the initial excitation of the disc. 
Planet formation models suggest that planetesimals 
should initially have very low eccentricities and incli- 
nations, damped by the gas disc. For this reason we 
choose to restrict our initial eccentricities and incli- 
nations to e max = 0.1 and I max = 10° in this work 
and allow any excitation above these values to occur 
naturally during the course of the simulations by in- 
teractions with the planet(s). This is a already a rel- 
atively _JrigJiJimi^£ojri£ar^ tomany observed debris 
discs dThebault & Augereaull2007t) . However, this ig- 
nores any potential excitation of the eccentricities and 
inclinations of the particles during terrestrial planet for- 
mation, after the gas disc has dissipated. A different 
choice of e max and I max could significantly alter our re- 
sults. Fig. [8] shows that the rate at which material is 
scattered inside of r zoc u increases with e max and Imax- 
To calculate these results it was necessary to start the 
simulations with particles spread initially over the full 
belt, so as not the bias the results, given that there is 
an increase in the size of the chaotic zone with both 
the initial eccentricity and inclination of the part icles 
dMustill & Wvattl 120111: IVeras & Armitagd |2004|) not 
accounted for in the formula of IWisdoml (ll980h ~This 
increase in the size of the chaotic zone, is responsible 
for an increase in the number of particles scattered and 
thus an increase in the rate at which particles are scat- 
tered inwards. 



4. Multi-planet Systems 

A more realistic planetary system is one in which there 
are multiple planets and the outer planetesimal belt is 
significantly further from the star. There is good evi- 
dence for the prevalence of multi-planet systems with 
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Label 



Mini 



1M 



Nep 



1M M , 
1M, 



2 



a p/ + 26a chaos 0.1 10° 
a p/ + 2Sa chaos 0.1 10° 
a p i + 2Sa chaos 0.1 10° 



/,„ or Using the same procedure as for the single planet 

7jgase (described in Sec. [3), the rate at which parti- 
IQgJes are scattered inside of r znc n = 1AU of the star in 
10 gur N-body simulations is calculated. Again, each par- 



e =02 
&max 0.3 CI 



c max 



dpi 
dpi 
Pi 



0.4 



pi 



<t p l 



,2 
-'")./ 

_2_ 



10- 4 

0.2 
0.3 
0.4 



10 



Full 



dpi 
a p i + 
dpi 



a p i + 2Sa c 



0.1 
0.1 
0.1 



S'cte is assigned a fraction of the outer belt mass and 
le rate at which material is scattered inwards is cal- 
Elated by considering the time between adjacent scat- 
lOPiring events. Again, we arbitrarily consider outer belts 
40° lOfKht extend to r be i, iOU , = \r be i uin = 93AU. We vary the 
10° 94glanet masses, M p \ and separation ARh between our 
10° 87^mulations. The dependence on the stellar mass M* and 
10° lOtte position of the outer belt, r\, e u,in should also be noted 
ttrrscussed further in Sec. |5j. The results are shown in 



0° 
20= 
30° 
40° 



Table 2. Details and results of the various one planet 
simulations. All simulations are for a \M^ ep mass 
planet, unless otherwise stated. 



a wealth of different architecture^. The full parame- 
ter space is clearly too large to investigate in this work, 
even if we restrict ourselves to N p i planets on circular 
orbits, in the same plane, interior to a planetesimal belt. 
For simplicity, we therefore focus on the example of a 
chain of equal mass planets, separated equally in units 
of their mutual Hills radii. This system should be con- 
sidered a representative system, from which it is possi- 
ble to better understand the scattering processes and the 
way that they apply to other systems. 

The simulations are set-up in the same manner as 
described in Sec. [2] We place the outer belt at rb e i r j„ = 
62AU, to match with the observations o f the outer belt 
in the Vega system dMiiller et al.l l2010h . although we 
anticipate that our simulations can be scaled to a lim- 
ited range of different belt radii. In the same manner as 
for the single planet system, the outer planet is placed 
such that it sculpts the inner edge of the outer belt, at 
dpi,out — fbeitjn - 5a c i mos . Equal mass planets are then 
placed interior to this, separated by ARh- In order to 
ensure that the planetary systems considered actually 
scatter material inside of r :oi a = 1AU, we consider the 
conservation of the Tisserand parameter, and the limit 
this provides on how far particles can be scattered by 
a planet (Eq.|6]). This time, we anticipate that the scat- 
tering by multiple outer planets excites the particles ec- 
centricities and inclinations before they are scattered by 
the interior planet. Therefore, we place a limit on the 
interior planet such that all particles with e > 0.2 can 
be scattered inside of r zo rf,=lAU. This implies that the 
planet should be at 5AU. Therefore we 'fill' our plane- 
tary system with planets equally spaced in terms of their 
Hill's radii down to a /,// ,„=5AU. The number of planets 
required to 'fill' the planetary system ranges between 3 
and 10 (as can be seen in Table©. A choice of a larger 
value for afmj„ could reduce the number of particles 
scattered inside of r zol a to zero. It is therefore worth not- 
ing that our simulations are more representative of the 
maximum, than minimum scattering rates. 



3 162/580 multi-planet systems according to exoplanets.org 
on 18 June 2012 



Fig- EI f° r the different simulations listed in Table|4] 

Fig- Hjclearly shows that all simulations follow the 
same general trend; the scattering rates decrease with 
time. Above this, there are two clear dependencies on 
the architectures of the planetary system, one on the 
planet mass and one on their separation. Firstly, as 
discussed in Sec. 13.21 we expect higher mass planets 
to scatter more material on shorter timescales. This is 
clearly seen in Fig. [9] where the scattering rates increase 
with planet mass between the purple blue crosses and 
yellow crosses. Secondly, the more widely spaced the 
planetary system, the lower the scattering rates e.g. the 
decrease in scattering rates between the yellow crosses 
(&Rh) and red crosses (15Rh). This follows simply 
from the fact that the wider the separation of the plan- 
ets, the smaller the proportion of the particles scattered 
by one planet that are on orbits that interact with the 
next inner planet, and so forth for each pair of planets. 
The planet mass clearly dominates the scattering rates. 

4.1. Application of these results to any system 

The simulations presented in this section are for very 
specific planetary systems, with strict constraints on its 
architecture. This was necessary in order for the current 
investigation to be feasible. Clearly, these do not mimic 
many, if any, real planetary systems. In fact, the current 
observations of multi-planet systems find a huge variety 
of planet masses and separations. Nonetheless, the re- 
sults from these simulations can be used to make useful 
conclusions about the scattering rates in more realistic 
planetary systems. Firstly, having understood the broad 
dependences of the scattering rates on planet separa- 
tion and mass, we are able to conclude that the highest 
scattering rates, at late times, occur for tightly packed, 
low-mass planets. Thus, our simulations with chains of 
Neptune mass planets separated by 87?// are representa- 
tive of some of the highest possible scattering rates af- 
ter millions of years of evolution, beaten only by chains 
of lower mass planets, packed even closer to the stabil- 
ity limit. We therefore consider that our simulations are 
broadly representative of an approximate upper limit on 
the scattering rates possible in any planetary system, al- 
though, it is worth noting, as will be discussed in Sec.|5J 
that we may have missed a few systems with very spe- 
cific architectures that enhance the rates still further. 
Secondly, the results presented here can be scaled, such 
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Planet 
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Table 3. A summary of the parameters for the chains of equal mass planets used in the simulations of Sec. 21 
shown in Fig. [9] The bottom half of the table describes scaled versions of these simulations that could be applied 
to Vega or rj Corvi. 
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Fig. 9. The rate at which particles are scattered, by a 
chain of equal mass planets of \M^ ep (dark purple) 
lMsat (blue) or IMj (red/orange), spaced equally in 
units of their mutual Hill's radii by either 8 (crosses) 
or 15 (stars) Rh- The rate is plotted as a fraction of 
the outer belt mass, in a belt spanning from 62AU to 
93 AU, and particles are scattered inside of r- 0£ /, = 1 AU, 
as described in Sec. |4] and Table |4] These results can 
be scaled to different systems, as described in Sec. 14.11 
The scattering rate is higher on shorter timescales with 
higher mass planets that are closely spaced, whilst on 
longer timescales tightly packed low mass planets scat- 
ter at a higher rate. 



that the scattering rates in a system of a known archi- 
tecture can be estimated. In this manner the scattering 
rates in a real planetary system, with planets on known 
orbits, could be estimated. 



In this work, we investigated the change in the 
scattering rates between chains of equal mass planets, 
evenly spaced in units of their Hill's radii (ARh), with 
different planet masses (M p i) and separations (A). To do 
this we fixed the stellar mass (M»), the radial location 
of the outer belt, (rt, e it ,,-„, rb e it,out), the radial location of 
the exozodi {r zoc a) and a parameter that we term a 
that specifies the number of planets, by providing an in- 
ner radial limit on their possible orbits. Real planetary 
systems are unlikely to have the values used for these 
parameters. We therefore discuss the effects of chang- 
ing the different parameters on the scattering rate, in 
such a way that we can understand and anticipate the 
behaviour in any planetary system. 

Firstly considering planetary systems of the form 
described in this section, with equally separated, equal 
mass planets, the scattering rates that we calculated (i.e. 
shown in Fig. [5J, can be applied almost exactly to an- 
other system in which the orbital timescales of the plan- 
ets/planetesimals are different. In other words, systems 
where the radial location of the outer belt or the stellar 
mass is changed. In this manner the results of our sim- 
ulations can be applied to planetary systems with outer 
planetesimal belts at any location, with the assumption 
that the radial location of the exozodi scales with the 
radial location of the outer belt. If the radial location 
of the outer belt is scaled by a factor, /, then the same 
scattering rates can be applied to the system by scaling 

3/2 3/2 

the time axis by a factor of aJ lout ~ r bell in . Similarly, for 
stars of different stellar mass, there is a small change in 
the scattering timescales that can be found by scaling 

1 /2 

the time axis by a factor M„ ' . 

There are three other parameters that can affect the 
scattering rates; the width of the outer belt (r/, e /, „,), the 
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location of the exozodi (r w dd an d the innermost planet's 
semi-major axis, a/m. The width of the outer belt is in 
essence unimportant, as long as it is greater than the 
planet's chaotic zone. It will however change the values 
for both the total outer belt mass, and the fraction of it, 
that is scattered inwards, in a manner that can be readily 
calculated using Eq.[3] 

The location of the exozodi (r zoc u) is crucial. 
Scattering rates generally scale with r :oc n, unless r zo< n is 
too small that no (or very few) particles can be scattered 
sufficiently far inwards (see discussion of the analytical 
limit in Sec. 13. U . Changing the location of the exozodi 
produces differences in the scattering rate of the same 
order of magnitude as the difference caused by chang- 
ing the planet separation (A). 

On the other hand, the value of the parameter a fm 
used here, was chosen to ensure that particles are scat- 
tered to the location of the exozodi, r :oi n- In a real plan- 
etary system, the innermost planet could take any po- 
sition, and the scattering rates will depend critically on 
this position. Test simulations find that scattering rates 
do not vary significantly with a fm, or equivalently the 
position of the innermost planet, so long as the planet 
separation does not increase. This only applies if a fm is 
not decreased so much that the inner planet is within 
a Hill's radius of the exozodi, or increased so much 
that particles can no longer be scattered from a/m to 
f Z odi- For the example system used in this work with 
dbeit,out = 62r zo di, analysis of our simulations found that 
particles are generally excited up to eccentricities of 
around e ~ 0.3 before being scattering by the planet. 
Thus, using Fig. [5] we anticipate that increasing a/m 
up to 10 AU, does not significantly alter the scatter- 
ing rates. We verified this with test simulations. In other 
words, taking the example of the chain of 4 Jupiter mass 
planets separated by 8Rh, the innermost planet is not re- 
quired in order to retain the scattering rates at the calcu- 
lated values. However, if a/m were to be increased be- 
yond 15AU, no particles would be scattered inwards. It 
is therefore easy to design a planetary system in which 
the scattering rates tend to zero, by applying Fig [5] to 
that system. 

Now, we consider the broad outcomes of these re- 
sults and their applicability to any planetary system. 
In general, the results presented in this work only fur- 
ther our understanding of planetary systems with a sin- 
gle, outer planetesimal belt and planets on circular, co- 
planar orbits, interior to the belt. Let us first consider 
chains of planets of non-equal mass. The outer planet 
controls the amount, and the rate, at which material 
is scattered into the inner planetary system. Thus, the 
outer planet dominates the scattering rates, whilst the 
inner planets modify these rates slightly by controlling 
how long it takes for particles to be passed inwards 
to the exozodi. This timescale depends on the mass of 
the planets. The best way to illustrate this is to con- 
sider the simple example of a chain of Neptune mass 
planets, that ordinarily scatter a small amount of ma- 
terial inside of r zoc n on long timescales (see Fig. [9]). If 
all but the outer planet were replaced with higher mass 
{e.g. ~ IMj) planets, whilst retaining the same separa- 



tion in Hill's radii, then material would continue to be 
scattered out of the outer belt on the same timescales, 
but would be passed from the outer belt to r l0 n more 
swiftly. Therefore, the rate of scattering of material in- 
side of r z „di would fall off more steeply with time. If the 
interior planets were of lower mass than the outer plan- 
ets, they could do the inverse and cause a bottleneck, 
delaying the scattering process. 

Now, all that remains in order to be able to apply 
our results to any planetary system, is a consideration 
of the effect of changing the separation of the plan- 
ets. Although our simulations clearly show that more 
widely spaced systems scatter at somewhat lower rates, 
as discussed in Sec. [4] the effect on the scattering rates 
is not as significant as that caused by changing the 
planet mass. The separation of the planets can, how- 
ever, be increased so much, that a gap in the planetary 
system is produced, over which very few particles are 
scattered. The separation at which this occurs can be es- 
timated by considering when the planets are so widely 
separated that q , nm(Tj) < a p i - AR H , or referring to 
lBonsor& WvattN2012[) for a more detailed calculation. 
To take an example, for a chain of Jupiter mass plan- 
ets, if any pair of planets were separated by more than 
20Rh, very few particles could be passed between them. 

To summarise, the arguments presented here can be 
used to calculate exactly the scattering rates in plane- 
tary systems that follow the mould of our example and 
to estimate the broad behaviour in any other planetary 
system. 

4.2. The collisional evolution of the outer 
planetesimal belt 

The outer planetesimal belt is a collisionally evolving 
system, like any debris disc. We have not taken the col- 
lisional evolution of the disc mass into account in our 
calculation of the rate at which material is scattered in- 
wards, shown in, for example, Fig. [9] Each individual 
data point in this plot was calculated by assigning each 
test particle a fraction of the total outer belt mass at the 
instant it was scattered inwards. Of course, if the disc is 
collisionally evolving, the total outer belt mass could 
potentially be different for each particle. Taking this 
into account, the new version of Fig. [9] would have an 
even steeper fall off with time. For practical purposes, 
the data in Fig. [9] will generally be used to assess the 
total outer belt mass required in order to obtain an ob- 
served scattering rate. This should be close to the disc's 
current mass, because in the general the timescales for 
the collisional evolution are not significantly faster than 
the scattering timescales. 

If we make the assumption that material must be 
scattered inwards at a rate of at least 10~ 9 M ffi /yr to re- 
supply an exozodi, the simulations shown in Fig.[9]can 
be used to determine the minimum mass required an the 
outer belt at 62 AU. It is important to consider whether 
or not such a mass can survive against collisions, at the 
given radial location, for the age of the system. This 
can be done analytically, for example using the mod- 
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els of IWvatt et all (l2007al) and Eq. [5] In fact, given that 
the results of our simulations can be readily scaled to 
systems with outer belts at different radial locations (as 
discussed in Sec. 14. Il l, the minimum mass required can 
be determined for belts at a variety of different radii. 
This scaling necessarily assumes that the location of 
the exozodi scales with the outer belt. The best way to 
illustrate this is to take the scattering rates, R s , calcu- 
lated and shown in Fig. [9] The disc mass required to 
produce these scattering rates, by an outer belt at any 
(reasonable) radial location, rb e it,m (i- e - scaling the rates 
of Fig. [9] with the orbital timescales), is given by 

RJM 9 yr~ l \ I r behMAU \ 3 / 2 
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We now have a requirement on the minimum mass 
needed in an outer belt, at a range of radial locations, if 
the outer belt is to supply material scattered inwards at a 
rate of at least l(T 9 M e /yr. Contrasting this with the col- 
lisional evolution of the disc, it is possible to determine 
which of these discs can survive against collisions for 
the age of the system. In other words a minimum limit 
on the radius required for the outer belt can be placed. 
This limit is shown in Fig. [10] Given the assumptions, 
this means that an outer belt that resupplies an exozodi 
must lie beyond the radial limit shown in Fig. [10] Each 
cross in this plot is taken from the simulations shown 
in Fig. [9] where the minimum mass required in the disc 
is calculated for a range of radial locations using Eq. [7J 
and the minimum possible radial location determined 
by substituting these masses into Eq.[6] to obtain: 
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Fig. 10. The minimum radius of an outer belt, rb e it,min, 
in which sufficient mass to scatter material inwards at 
a rate of 10~ 9 M e /yr can survive against collisions, as 
a function of time. Each cross is taken from the simu- 
lations shown in Fig. [9] scaled to different radii for the 
outer belt and planets using Eq. [7] and converted to a 
minimum radial location for the belt using Eq.[9] 



4.3. The detectability of the outer belt 

Our simulations can be used to place limits on the re- 
quired masses and radial locations of an outer planetes- 
imal belt, in order that it could potentially scatter mate- 
rial inwards at sufficiently high rates to retain the cur- 
rently observed levels of dust in exozodi. These limits 
can be compared with our ability to detect such discs. 
Taking our simulations for chains of multiple equal 
ass planets, shown in Fig. [9] we calculate the mini- 
um disc mass required to scatter material inwards at a 
fiducial rate of 10~ 9 M ffi /yr, at an age of lOOMyr. This is 
the total mass, in bodies of all sizes, in a disc between 



(0^/^)3/1300/^-3)3/13(^)3/13(^/^-1)5/26.^^ + ^ ^ 



3r,». 



This makes the assumption that the belt width is 
small and takes the values for the parameters listed in 
Sec. 12.31 Fig. [TOl clearly shows that because relatively 
massive planetesimal belts are required to scatter ma- 
terial inwards at the required rates, the belt must be 
reasonably far from the star in order that such a mass 
can survive against collisions for the typical ages of ex- 
ozodi systems. This of course assumes that the colli- 
sional evolution of the belt commenced very early in 
the star's evolution and that all bodies in the disc are in- 
volved in this collisional evolution. It also assumes that 
the belt width is small and fixed at — = 0.5, that the 

r 

location of the exozodi scales with the radial location 
of the outer belt and that a scattering rate of 10~ 9 M ffi /yr 
is reasonable requirement to retain the observed dust 
levels in an exozodi. In this case, Fig. [10] shows that 
in order for an outer planetesimal belt to scatter dust 
inwards at this level for longer than lOOMyr, it must be 
outside of ~ 20AU, otherwise insufficient mass can sur- 
vive against collisions. This is an important conclusion 
because many observed debris discs are inside of this 
radius. 



As discussed 

in Sec. 14. II the scattering rates, and thus the minimum 
required disc mass, can be scaled to outer belts at any 
radial location, although only with the assumption that 
the position of the exozodi, r zoc n scales in a similar man- 
ner. Thus, every point in Fig.[9]can be considered to rep- 
resent the scattering rate after lOOMyr in a belt at a dif- 
ferent radial location. The radial location and minimum 
required masses for the outer belts to scatter material 
inwards at rates higher than 10~ 9 M e /yr are shown by 
the blue crosses in Fig. Q~T] for a chain of \Ms a t plan- 
ets separated by 8Rh- This is the simulation with the 
largest detectable parameter space of any of our simu- 
lations. This plot reiterates the fact that massive, large 
radii planetesimal belts are the most adept at scattering 
material inwards. 

Some constraints on the existence of the discs out- 
lined on Fig.QT]can be made by considering their col- 
lisional evolution. As discussed in Sec. 12.31 only belts 
below a maximum mass, can survive at a given radial 
location for a given age. This maximum mass is plotted 
by the dotted line on Fig. [TT] calculated usi ng Eq. [5] 
taking the parameters of IWvatt et al.l d2007bl) . listed in 
Sec. 1231 
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The next question to consider is which of these discs 
could be detectable. One way to detect a debris disc is 
by an infrared excess, for example at 70ji/m. For this 
plot, we consider Spitzer's detections limits, as Spitzer 
surveys of debris discs are comprehensive and well 
understood, despite the better detection capabilities of 
Herschel or ALMA (for colder debris discs). Spitzer 
observations are calibration limited, requiring a ratio of 
the disc to stellar flux above > 0.55 for a detection. 
By considering a simple model for a debris disc, with a 
steady-state collisional cascade, n{D)dD oc D 1 ^dD, for 
grai ns between th e blow-out size, Dm and D BMl = 60km 
(e.g. lWvatd2008h . the minimum detectable disc flux can 
be turned into a minimum detectable disc mass, Ma et , 
of: 



M de , = 1.6 x 10 10 J? M ,« m ^(D b ,/fi,n)(D n 



Jkm) 



(9) 



(_LJLe_\ B V (A,1 
\(T,/Kf) B v O,T, 



'disc) 



where the B v refers to the Wien function describing 
black-body emission at a given wavelength, A or fre- 
quency v, at either the stellar temperature, T* or the 
disc temperature, Tn sc . Discs that can be detected are 
those than lie above the dashed line in Fig.Qj] for a G2 
(1M©) star. The dark grey shaded area shows all discs 
that could be detected by Spitzer and scatter material 
in at a sufficiently high rate to resupply an exozodi at 
10~ 9 M ffi /yr, at an age of lOOMyr. The lighter shaded 
area shows the parameter space in which debris discs 
capable of scattering material inwards at a sufficiently 
high rate could lie, but not be detected with Spitzer at 
70jum. 

The clearest conclusion from this plot is that debris 
discs exist that could supply scattering at the required 
rates in systems older than lOOMyr. Such debris discs 
are at large radial distances and massive. It is evident 
from Fig. [TT] that such large radii discs are the popula- 
tion about which least information is known, although 
this may change in the near future with instruments 
such as ALMA or SCUBA 2. 

This plot should be viewed with several uncertain- 
ties and restrictions in mind. Firstly, all the lines on 
the plot vary with the stellar luminosity or mass of 
the star. This plot assumes that required scattering rate 
is 10~ 9 M ffi /yr, a value which is uncertain even for the 
Vega system and may vary significantly between dif- 
ferent systems. It also makes the potentially unrealistic 
assumption that the radial location of the exozodi scales 
with that of the outer belt. 



4.4. Application to observed systems 

The number of systems with detected exozodi is lim- 
ited, only a handful of the detected systems also have 
outer planetesimal belts and there are very few planet 
detections. For these reasons, we decide to focus here 
on two example systems, both of which have resolved 
outer belts, and detailed modelling of the inner exozodi. 
Neither of the systems have any evidence for planetary 



1000.00 
100.00 
10.00 
1.00 

0.10 
0.01 



Max / 

/ 

/ 

/ 

/ G2 



1 



10 100 

Radius / AU 



1000 



Fig. 11. Outer belts capable of scattering material in- 
wards at 10^ 9 M e /yr after lOOMyr, around a G2 (1M Q ) 
star, with a chain of \Ms a t planets separated by %Rh, 
that can be detected by Spitzer at 70/vm are shown in 
dark grey, whilst those belts that cannot be detected are 
shown in light grey. This is the lowest limit of the simu- 
lations shown in Fig. [4] The dotted line shows the maxi- 
mum mass that can be retained against collisions (Eq.[5j 
for lOOMyr and the dashed line shows the Spitzer de- 
tection limit for a G2 star. Full details of this plot are 
described in Sec. 14.31 



companions. The two example systems were chosen to 
contrast. Vega is a 450Myr old A star, whilst 77 Corvi 
is an older (Gyr), sun-like star, with an extremely large 
infrared excess. We apply the simulations presented in 
this work, to these systems, in order to determine the 
possible scattering rates that could occur, if we make 
the hypothesis that chains of planets orbit between the 
inner and outer dust discs. We then compare these scat- 
tering rates with observations of the exozodi and outer 
belt, and assess the feasibility of such scattering to solve 
the conundrum of the high levels of exozodiacal dust 
observed. 



4.4.1. Vega 

Vega's outer belt was the fi rst debris disc to be de- 
tected ( lAumann et al.l 1 19841) and has since been ob- 
served many times at a variety of different wavelengths. 
The exact properties derived from observations of the 
outer belt vary in t he literature, but he re we consider 
the best fit model o f lMuller etalJd20Toh . that takes into 
account observations from the near-infrared through 
to sub-mm wavelengths. Their best fit model finds a 
dust belt between 62 and 120 AU, with 5.09M ffi of 
dust and a total mass of 46.7M e , which collisionally 
evolved from an initial mass of 55.5M e . Fitting of 
isochrones finds a n age of 455 ±13 Myr and a mass of 
2.157 ± 0.017M e dYoon et alJIiOlOh. Exoz o diacal dust 



was dete cted around Vega by Absil et ail (120061) and 
again by IDefrere et all d201 lb . The best-fit to the ob- 
servations finds a mass of 10~ 9 M m in small dust grain s 
between 0.2/im and 1mm in size (IDefrere et al.ll20lTb . 
A collisional timescale of ~ lyr is derived. This sug- 
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gests that if the disc is to be retained at its current levels 
it needs to be replenished at a rate of at least 10~ 9 M e /yr. 

The scattering rates determined in our multi-planet 
simulations described in Sec. |4] and shown in Fig. [9] 
can be applied to Vega. Without any limits on plan- 
etary companions less than several Jupiter mass (e.g. 
iHinklev et aT]|2007t lltoh et al.ll2006l) . we are free to hy- 
pothesise the existence of any inner planetary system. If 
Vega is really orbited by as yet undetected planets, they 
may or may not resemble the chains of multi -planet sys- 
tems considered in this work. However, by considering 
the hypothesis that there are chains of equal mass plan- 
ets, interior to Vega's outer planetesimal belt, we are 
able to place an estimate on the maximum rates that 
planets could possibly scatter material inwards. 

The position of the outer belt was already chosen to 
fit with that of the outer belt in the Vega system. Thus, in 
order to apply the results of these simulations to Vega, it 
is only necessary to scale the scattering times by a fac- 
tor of , 1 ~ 0.68 in order to account for the increased 

V2.157 

stellar mass. These results now represent systems with 
chains of equal mass planets still separated by A = 8/?// 
or 15Rh, but with increased masses of IM^ep, 2M sa , or 
2Mj. This is detailed in Table|4] In Fig.Q~2]we calculate 
the total mass in the outer belt if these chains of equal 
mass planets were to scatter material in to the position 
of the exozodi, which we take to be 1 AU, at the required 
rate of 10~ 9 M e /yr. This is the minimum required mass 
in the outer belt, assuming a constant surface density of 
I.(r)dr oc r~ l dr, between 62 and 120AU. 

Fig. Q~2] clearly shows that by the age of Vega 
(455 ±13 Myr) a very massive outer belt is required 
in order to sustain the high levels of dust observed. In 
general the required masses are higher than the obser- 
vationally derived mass for the belt (shown by the dia- 
mond), or ev en the pre-collision al evolution initial mass 
of 55.5M ffi (iMuller et al.1 l2010h . Given that the plot- 
ted values assume that all material scattered inwards is 
converted to the observed small dust grains and thus 
are likely to underestimate the belt mass required, this 
places doubt on the ability of Vega's outer belt to re- 
plenish the exozodi. Such a conclusion, should, how- 
ever, be viewed in the light of the large uncertainties, 
in particular on the rate at which Vega's exozodi needs 
to be replenished. A small decrease from the currently 
used value of 10~ 9 M e /yr, would mean that a chain of 
low mass planets could readily scatter material at the 
required rates. Another possibility is that the scattering 
process in Vega started less than 455Myr ago and there- 
fore that the scattering rates are higher at the current 
epoch. To summarise, the uncertainties on the values of 
various parameters are sufficiently large that these sim- 
ulations by no means rule out steady-state scattering by 
a chain of planets to replenish Vega's exozodi, however, 
the evidence presented in favour of such scattering re- 
quires the existence of many closely spaced, terrestrial 
to Saturn mass planets, orbiting between the two belts. 
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Fig. 12. The total outer belt mass required if scattering 
by a chain of equal mass planets, as shown in Fig.|9l or 
detailed in Table [4] is to replenish an exozodi inside of 
1AU around Vega at the required rate of 10~ 9 M ffi /yr. 
Each cross represents an individual scattering event 

shown in Fig. [9] and scaled by a factor of ( 2 '\ 5 ™ 9 ) 
to account for Vega's stellar mass. Vega is plotted by 
the d iamond at its age of 455 + 13 Myr dYoon et al.1 
l2010h and with an o uter belt mas s of 46 .7M ffl , as found 
by the modelling of IMuller et"al] d2010l) . No error bars 
are quoted on the mass, however, its model dependence 
should be noted. 



4.4.2. 77 Corvi 

7] Corvi as a significantly older (1.4 + 0.3Gyr) system, 
with an extremely high infrared excess. Observations 
of 77 Corvi in the mid-infrared and sub-mm sug- 
gest that it has both a cold and a warm compo- 
nent dSmith et al.ll2009bt IChen et al.l|2006h IWvatt et al l 



120051) . Both lSrrrith et alJ d2009bl) and fLisse et all (I2012T 
find that th e warm component is inside of 3.5AU and 
iLisse et ail d2012l) find that it can be explained by 
9 x 10 18 kg ~ 10 _6 M e of dust in particles between 0.1 
and 100/jm in size, with a lifetime of 3-6 yrs against 
collisions. This suggests that this dust needs to be re- 
plenished at a rate of ~ 10~ 6 M e /yr in order to retain 
the current observed leve ls. Sub-mm obser vations of 
the outer belt, modelled in lWvattetai]d2005l) find a ra- 
dius of 150 ± 20AU and a mass of ~ 20M e . 

Again the scattering rates shown in Fig.|9]can be ap- 
plied to rj Corvi, hypothesising the existence of chains 
of equal mass planets orbiting interior to the outer belt. 
Again, the planetary system of 77 Corvi may differ radi- 
cally from such chains of equal mass planets, nonethe- 
less, such systems provide a useful upper limit on the 
possible scattering rates. In this case, the scattering 
rates determined in Fig. [9] apply at times greater by a 

factor ( 1 §m-) V2 (rm;) V2 ~ 2.57, to account for the 
differences in belt radius and stellar mass. The scatter- 
ing rates now apply to chains of planets with slightly 
higher masses, scaled upwards by a factor of 1 .4, as de- 
tailed in Table |4] The scattering rates now apply to a 
system with an exozodi at 2AU, slightly inside the ob- 
served position of less than 3.5AU, however, this could 
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Fig. 13. The same as Fig. [12] but for r\ Corvi, with its re- 
quired scattering rate of ~ lCT 6 M ffi /yr. The scaling fac- 
tor this time is (^^) 3/2 (r^) 1/2 , to account for the 
difference in stellar mass and radius of the outer belt. r\ 
Corvi is plotted by the cross, at its age of 1 .4 ± 0.3Gyr 
and w ith an outer belt mass of ~ 20M m ( IWvatt et al.l 
2005]). A gain, no error bars are quoted on the mass 
of the outer belt, but its model dependence should be 
noted. 



only potentially cause a small change in the calculated 
scattering rates. 

The level of excess observed in this system, and 
therefore the replenishment rate required, is signifi- 
cantly higher than for Vega's exozodi. This means that 
the minimum required outer belt mass to replenish the 
exozodi at ~ 10~ 6 M ffi /yr, shown in Fig. [13] are ridicu- 
lously high. This is clear evidence that the warm dust 
observed in 77 Corvi is at levels far too high for it to be 
resupplied by the scattering of material inwards from 
the observed outer planetesimal belt by a chain of stable 
planets. Another explanation is required for this system. 



5. Discussion 

In this work we have considered some example plane- 
tary systems in order to investigate the hypothesis that 
high levels of exozodiacal dust observed in some plan- 
etary systems could be retained by material scattered 
inwards from an outer planetesimal belt. We used N- 
body simulations to determine the scattering rates from 
chains of equal mass planets, on circular, co-planar or- 
bits, interior to an outer planetesimal belt. These were 
used to assess the validity of this hypothesis by com- 
parison with observations. We used our simulations to 
outline the properties (mass and radius) of debris discs 
that are capable of scattering material inwards at suffi- 
ciently high rates, which were then compared with our 
ability to detect such discs and whether the masses are 
consistent with collisional evolution. In general, this is 
highly dependent on the properties of individual sys- 
tems and thus, the most definitive conclusions can only 
be made by comparison with observed systems, as was 
done here for Vega and 77 Corvi. We now discuss some 



of the uncertainties that should be taken into account 
when considering the results presented in this work. 

5.1. Uncertainties 

The biggest uncertainty in these comparisons is the de- 
termination of the rate at which the exozodi must be 
supplied with material in order to retain the observed 
levels of exozodiacal dust. This is calculated from a 
fairly uncertain mass of dust in the exozodi, combined 
with an even more uncertain timescale on which the 
dust is removed. The mass calculated for the dust is 
very model dependent, with a particular dependence on 
the size distribution and s ize range of g rains consid- 
ered. For Vega's exozodi, iDefrere et al.l d201 ll) create 
a grid of model discs, varying the composition, optical 
properties, radial and size distribution of the dust, and 
then use Bayesian statistics to assess which model pro- 
vides the best fit to the observations. The largest grain 
size they consider is 1mm, but this upper limit is poorly 
constrained. The disc will be removed either by radia- 
tive forces, on dynamical timescales of ~yr, or on col- 
lisional timescales, also ~yr. Both these timescales de- 
pend on the exact properties of the disc and are fairly 
uncertain. The radial location of the dust is also not pre- 
cisely constrained, however, this does not significantly 
effect the scattering rates. 

Similar uncertainties, although less critical in terms 
of the current hypothesis, also apply to our knowl- 
edge of the properties of the outer disc. For the sys- 
tems where the outer belt is resolved, the radial lo- 
cation of the dust belt is fairly well constrained. 
The total mass of the disc, however, depends on 
the extrapolation of the observed population of small 
dust particles to larger sizes. Although model depen- 
dant, detailed numerical modelling (e.g.lThebault et all 
" Jjhebault & Augereaull2007t 



2003 



JKrivovet all 120061: 
Lohne et alj|2008l) finds that the ratio of the total mass 
of the disc to the observed dust mass calculated using 
realistic size distributions does not differ significantly 
from simplistic estimates made using the size distribu- 
tion from a steady state co llisional cascade (n(D)dD oc 
D 1/2 dD dDohnanvi) fT969h with a power index of 
?6) 



dTanaka et al. 1996b ). Thus, the calculated values for 



the mass should not vary too dramatically. The biggest 
uncertainty is in the size of the largest body in the colli- 
sional cascade, which can alter the calculated mass sig- 
nificantly. It is also possible that larger bodies exist that 
are not yet in collisional equilibrium. These could con- 
tribute to the scattering process, but only in a stochastic 
manner. Given our total lack of information regarding 
such bodies, it seems unreasonable to invoke their con- 
tribution to the scattering rates, although, evidence ex- 
ists in some systems that such stochastic collisions ma y 
produce the observed exozodi dLisse et a l. 2009, 2012). 

The second biggest uncertainty in our calculations 
regards the conversion of the bodies scattered inwards 
from the outer belt, with a size distribution reflecting 
that of the outer belt, to the observed small dust. In 
order to obtain maximum possible scattering rates, we 
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considered that this process is 100% efficient. There are 
three main processes by which larger bodies are con- 
verted to small dust, collisions, evaporation or sponta- 
neous disruption. Marboeuf et al, in prep investigated 
in detail the evaporation of comets as a potential mech- 
anism to produce exozodiacal dust. They find that it is 
reasonably efficient in high luminosity systems, with a 
maximum possible efficiency that reflects the dust:ice 
ratio in the comets (generally about 1, so 50% efficient). 
Dust may also be produced as larger bodies are ground 
down by collisions, or cometary bodies may sponta- 
neously disrupt, as suggested as an explanation for th e 
Solar System's zodiacal cloud dNesvornv et aljBOlOl) . 
None of these processes will be 100% efficient at con- 
verting large bodies to the observed small grains. In this 
work, we merely note that where outer belt masses have 
been calculated, these are very much minimum values, 
rather than attempting to make any conclusions regard- 
ing the conversion of the scattered bodies to the ob- 
served small dust. 

There are then several uncertainties based on the 
initial conditions used in our simulations, which were 
discussed in detail in Sec. 13.31 It is, however, worth re- 
iterating a dependence of our results, on the fact that 
planetesimals remain on orbits such that they can be 
scattered by the planet, at the end of planet formation. 
Should this prove not to be the case, whether due to 
planetary migration, or an inaccurate estimation of the 
parameter that we used to characterise this clearing, 
«,„,„, our scattering rates will only decrease. Again, re- 
iterating the fact that the outer planetesimal belt masses 
determined by our simulations are very much minimum 
values. We also assumed that all planets were fully 
formed at t = 0, which of course is not generally true, 
particularly for terrestrial mass planets. During terres- 
trial planet formation, the scattering processes are likely 
to be similar to those seen in our simulations, albeit con- 
trolled by the somewhat smaller planetary core. Thus, 
we do not anticipate that our results would be signifi- 
cantly altered if planet formation was only completed 
some time after t — 0. 

Of potential concern is how representative our ex- 
ample systems are of the full range of real planetary 
systems. Such well ordered planetary systems are un- 
likely to be found (often) in nature. In particular, it is 
unlikely that planet formation would have been equally 
efficient throughout any planetary system. They should, 
however, be representative of the range of scattering 
rates possible from the full range of real planetary sys- 
tems, as a function of time. This hypothesis is supported 
by test simulations of two further planetary systems, 
with random architectures, which were found to scatter 
material inwards at rates that lie neatly within the range 
of scattering rates found for our example systems. As 
discussed, for example in Sec. 14. II or Sec.|4] it is easy 
to design planetary systems with 'blockages' where the 
scattering rates tend to zero. However, given our un- 
derstanding of the dependence of the scattering rates 
on the planet masses and separations (see discussion in 
Sec. 14. II ), we anticipate that very few systems can scat- 
ter material inwards after millions of years of evolution, 



at rates higher than those outlined here for a chain of 
tightly packed Neptune mass planets (given constraints 
on the outer belt mass, stellar mass, age etc). It is, how- 
ever, still possible that a few very specific architectures 
exist where the scattering rates are enhanced beyond 
those seen in our simulations, for example due to res- 
onant interactions. Alternatively systems with eccentric 
planets or planetesimals may prove to scatter material 
inwards at higher rates, although test simulations sug- 
gest that they can only do so on short timescales. 

There are two other potential contributors to the 
scattering rates, that we have ignored so far. Both con- 
tribute at such low levels that we feel justified in leav- 
ing them out of our calculated rates. The first is par- 
ticles that started on orbits between the planets. If the 
planets are sufficiently widely separated, it may be that 
some material was left orbiting between the planets, 
rather than being accreted onto the planets. In general 
such planetesimals are scattered on short timescales, 
however, if the separation of the planets is larger than 
their dynamical influence, i.e. ARh > 5a c i mos , planetes- 
imals may remain on stable orbits between the planets, 
even on long timescales. The scattering process for such 
planetesimals is very similar to that for the outer planet 
in our simulations, except on shorter timescales, given 
the shorter orbital periods. Thus, at late times the contri- 
bution of such particles is likely to be small. The only 
question therefore is whether such belts of planetesi- 
mals, or indeed planetesimals in the process of being 
scattered inwards, may be detectable. 

The second potential contribution is from par- 
ticles that started far outside of the outer planet's 
chaotic zone and are scattered on long timescales. In 
a multi-planet system, the dynamics are more com- 
plicated and there are a range of secular and multi- 
body resonances that can influence particles in the 
belt. These have been shown to provide a small flux 
of material i nto Neptune's chaot i c zone in our Solar 



System (e.g. Duncan etal.1 19951: iHolman & Wisdom 



1993 



Levison & Dunca ni l 19971: lEmel'vanenko et al 



2004t iMorbidellil 1 1 9971) . This contribution cannot be 



significant as in our simulations to test this effect with 
a particle density of 50 particles per AU, not a single 
particle, that started outside of a ma x = a p \ + 25a c i mos and 
inside of 100AU, was scattered inwards, although there 
was evidence of the effect of the inner planets on the 
outer regions of the belt. 



5.2. Comparison with real planetary systems 

All of the potential uncertainties discussed in the pre- 
ceding section need to be taken into account when com- 
paring the results of our simulations with observed ex- 
ozodi. In Fig. QT| we determined the mass and radius 
required for the outer belt, such that it can scatter mate- 
rial to the location of the exozodi at a sufficiently high 
rate. These requirements could change significantly, for 
example, if our estimate on the rate at which mate- 
rial needs to be scattered inwards was wrong, or the 
conversion of the scattered large bodies to small dust 
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is inefficient. Changes of several orders of magnitude 
are already seen between the example systems consid- 
ered here. Given a diversity of planetary systems, this 
would naturally result in a diverse population of exo- 
zodi, that could help to explain the range of exozodi 
detections and non-detections found so far. These un- 
certainties render insignificant the possible differences 
in the scattering rates caused by the unrealistic assump- 
tion that the radial location of the exozodi scales with 
that of the outer belt. 

A general conclusion from Fig.QT|is that the longer 
lifetimes and scattering timescales of large radii outer 
belts are required to sustain high scattering rates on 
long timescales. Unfortunately the population of large 
radii debris discs is the least well constrained, although 
sub-mm instruments such as ALMA and SCUBA2 may 
shed light on this in the near future. For a more detailed 
comparison, we considered two example systems, Vega 
(Sec. 14.4. U and tj Corvi (Sec. 14.4.21 ). For 77 Corvi, the 
mass of observed outer belt is orders of magnitude too 
low to retain the exozodi at its currently observed lev- 
els and it therefore seems more likely that there is an- 
other explanation for this system. For Vega, the replen- 
ishment of the exozodi by a chain of planets seems un- 
likely, but given the uncertainties discussed above, it is 
not ruled out, although very specific constraints on the 
system would be required. 

Our N-body simulations necessarily only consider 
a very limited range of all possible planetary system 
architectures. We consider these to be a fair represen- 
tation of the range of scattering rates possible in any 
arbitrary planetary system, and discussed their appli- 
cation to any planetary system in detail in Sec. 14.11 
However, it still remains possible that we have missed 
some very specific architecture involving unequally 
spaced planets of different masses, potentially on ec- 
centric or inclined orbits, that is significantly more ef- 
ficient at scattering material, for example where sec- 
ular or mean motion resonances enhance the scatter- 
ing rates. However, within the limitations of reasonable 
computing timescales, a full investigation was not pos- 
sible in this work. 

5.3. Application to unstable systems 

The results presented in this work find high scatter- 
ing rates for the first thousands or millions of years 
of a system's evolution, but much lower rates in stable 
planetary systems after hundreds of millions of years. 
However, there is a potential way to get high, albeit 
temporary, scattering rates in such older mature sys- 
tems. This could happen if the empty chaotic region sur- 
rounding a planet is refilled with fresh material. A way 
to achieve this is if the planet is displaced from its cur- 
rent orbit by a dynamical instability and scattered into a 
planetesimal belt, or if it migrates into such a planetes- 
imal belt at late times. Such processes could restart the 
scattering processes and significantly increase the scat- 
tering rates. From our simulations, we can estimate that 
the scattering rates could remain high, even for thou- 



sands to millions of years after such an event. However, 
this is only an estimate because we did not consider any 
test particles closer to the planet than «„„■„ = a p \ + <5 ",°" 
and if the chaotic zone were refilled with planetesimals 
after an instability, it is likely that the planetesimal belt 
will be cold. Further detailed work is required to inves- 
tigate this hypothesis. 



6. Conclusions 

In this work we have considered the conundrum that 
high levels of exozodiacal dust are observed in some 
systems, despite it being clear that such small dust 
has very short lifetimes against collisions or radiative 
forces. We considered the hypothesis that the exozo- 
diacal dust could be replenished by material from an 
outer belt, scattered inwards by planets. We presented 
N-body simulations that calculate the efficiency of scat- 
tering by a chain of equal mass planets, in a way that 
can be applied to a wide range of different planetary 
systems, irrespective of whether or not exozodiacal dust 
has been detected. The outer planet scatters planetes- 
imals as it sculpts the inner edge of the planetesimal 
belt. Particles are cleared from the chaotic zone on rela- 
tively short timescales, but these particles may continue 
to be scattered multiple times by the planet, evolving at 
constant pericentre or apocentre, until they are scattered 
sufficiently far inwards that they are passed onto an in- 
terior planet or ejected. In this manner particles can be 
passed along a chain of planets and may only reach the 
location of the exozodi, some after hundreds of millions 
of years. 

Our simulations show that it is possible for this scat- 
tering to be at a sufficiently high rate that it could retain 
the currently observed levels of exozodiacal dust. This, 
however, places very stringent constraints on the archi- 
tecture of the planetary system, namely that it contain 
a chain of closely spaced, low-mass planets orbiting in- 
terior to a massive, large radii outer planetesimal belt. 
Without a complete census of the architecture of plane- 
tary systems and/or detailed statistics about the exozodi 
population, it is difficult to make definitive conclusions 
about the proportion of systems in which this scenario 
could apply. We considered in detail two systems with 
observed outer belts and exozodi, and we hypothesise 
that there may be as yet undetected planets orbiting be- 
tween the observed exozodi and outer belt. We find that 
for 77 Corvi, an old, sun-like star with very high levels 
of infrared excess, that the exozodi cannot be retained at 
its current levels, on Gyr timescales, with this process, 
whereas for Vega a very contrived architecture would 
be required for the planetary system to retain the levels 
of exozodiacal dust observed around Vega. 

In general in this work we show that material is scat- 
tered into the inner systems of many planetary systems, 
at rates that depend on the architecture of the planetary 
system, but are in general low, particularly after hun- 
dreds of millions of years. We put this forward as an 
explanation for the observations of high levels of exozo- 
diacal dust, but find that sufficient material is only scat- 
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tered inwards in systems with very contrived architec- 
tures. An alternative solution for, particularly the older 
planetary systems with exozodiacal dust, could be that 
we are observing these systems in the aftermath, albeit 
the long aftermath, of a dynamical instability that led to 
a planet being scattered into the outer belt. 
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